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Abstract 

Essential oils have antifungal and antipathogenic effects and therefore are tar- 
gets in plant pathology research for their potential uses as natural substitutes 
for inorganic plant pesticides. Beauveria bassiana, an entomopathogenic fun- 
gus, can endophytically colonize a vast number of plant species and trigger 
induced systemic resistance against plant pathogens. Spore germination is 
the most vulnerable in the fungal life cycle and is therefore a good candidate 
for monitoring the effect of essential oils on the growth of B. bassiana. Per- 
centage germination of fungal spores and length of germination tubes were 
recorded from experiments. A mathematical model that was able to capture 
the effects of cymene, an essential oil produced by Monarda, on the germi- 
nation and growth was developed. This is the first report of a model for the 
impact of essential oils on B. bassiana spore germination. 
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1. Introduction 

Monarda, a member of the family Lamiaceae native to the Great Smoky 
Mountains National Park, synthesizes highly volatile substances, called es- 
sential oils. Essential oils have an antifungal and antipathogenic effect and 
therefore are targets in plant pathology research because of their potential 
as natural substitutes for inorganic plant pesticides. 

The role of essential oils in control of Rhizoctonia damping-off in tomato 
plant with bioactive Monarda herbage was the focus of the study conducted 
at the University of Tennessee in Knoxville. Tomato seeds grown in a 
medium containing essential oils extracted from Monarda were protected 
against damping-off disease caused by fungus Rhizoctonia solani [1]. Beau- 
veria bassiana, a beneficial, entomopathogenic fungus, can endophytically 
colonize a vast number of plant species and can also trigger induced systemic 
resistance against plant pathogens [2]. The overall goal of this project was to 
determine if the combination of the treatment of tomato seeds with bioactive 
herbage and with B. bassiana enhances the resistance of the tomato plant 
against its fungal and bacterial pathogens. 

A spore germination of B. bassiana occurs after a spore dispersed from 
a colony lands on a moist surface and is activated by conditions that are 
optimal for fungal growth, such as heat, moisture and nutrients (Fig. 1). 
As water molecules from the environment enter the spore, the spore starts 
to swell; its size increases until it reaches a critical volume when a small 
protrusion starts to form on one side of the spore. This protrusion will 
develop into a branching germ tube that will develop into the fungal hyphae. 
This initial stage of fungal growth is the most vulnerable in the fungal life 
cycle and is therefore a good candidate for monitoring the effect of essential 
oils on the growth of B. bassiana. 

The overall aim of this study was to model the effect of cymene, an 
essential oil produced by Monarda, on the germination and growth of the 
beneficial fungus, B. bassiana. Our specific objectives were to measure the 
spore germination percentage and model the spore germination tube growth 
rate as a function of cymene concentration. 

2. Materials and Methods 

Combinations of cymene and ethanol were used to produce six different 
concentrations of cymene: 500, 50, 5, 0.5, 0.05, and 0.005 mM. Ethanol was 
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Figure 1: Fungal spore germination depicting the stages of spore germination formation 
and branching 

used as the control. A 200 /il capacity pipette was used to extract 128.4 
/il of ethanol and a 10 /il capacity pipette was used to extract 10 /il of 
undiluted cymene to a small, previously marked test tube in order to obtain 
a 500 niM cymene concentration. A test tube shaker was used to stir the 
compound until it was homogeneous. A pipette tip was changed prior to 
being used in a newly obtained dilution to prevent mixing of two different 
concentrations. In a new test tube, 10 /il of 500 mM cymene concentration 
was diluted with 100 /il of ethanol to acquire 50 mM cymene concentration. 
Similarly, 10 /il of each newly obtained cymene concentration was diluted 
by 100 /il of ethanol to acquire ten times smaller concentration until 0.005 
mM cymene concentration is obtained. Each time the mixture was stirred 
by avartex shaker and pipette tips were changed prior to handling another 
cymene concentration. 

The environment for the fungal growth was set up in a Petri dish by 
placing two microscope slides on the opposite sides of the Petri dish and half 
of the filter glass in front of the microscope slides (Fig. [2]). After 200 /il of 
water agar was applied by a pipette to both microscope slides to provide water 
that is essential for fungal spores to grow, 10 /il of an oil concentration was 
applied to the filter glass. A sterilized wire spatula was used to put spores into 
a test tube filled with water. Since spores tend to cling together in water, a 
small amount of detergent (Tween 20) was applied to the dilution to disunite 
spores. The dilution was sonicated to homogenize mixture. Finally, 10 /il 
of the dilution was applied on the top of a water agar on each microscope 
slide. The Petri dish was closed and placed on a wet sponge in order to 
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Figure 2: Environment for growing Beauveria bassiana 



keep the spores moist over the time period of the experiment. The procedure 
was repeated until all the cymene concentrations and the pure ethanol were 
applied to three Petri-dishes each. The date and cymene concentration were 
indicated on each Petri dish. The experiment was performed twice. 

At certain time points, measurements were done for each cymene con- 
centration and control. The microscope slides were placed under microscope 
with a camera attached at the top. Multiple slides of spores were captured 
under the same magnification for each oil concentration. Since the micro- 
scope had a calibrated reticle for measuring, a distinct spore, with or without 
a germination tube, would be chosen and its length would be read from the 
microscopic ruler. After the images were printed, the chosen spore was found 
and its length was indicated next to it. Images were separated by concen- 
tration and time point; each concentration had a different reference spore of 
known length. The first data set was obtained at 16, 23, and 40 hours. The 
second set of measurements was performed at 12, 16, and 23 hours. When 
images from all concentrations at each hour were obtained and printed, ger- 
mination tubes were measured by a ruler and entered to Microsoft Excel 
spreadsheet in millimeters. Data sets, time points, and concentrations were 
differentiated in the spreadsheet, and the length was converted from millime- 
ters to micrometers. The formula used was: 

L * M 
1mm = — - — fxm 

where L represents length of the spore, read from the microscope, M repre- 
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sents the microscope magnification, and I represents length of the same spore 
in millimeters, measured from the print out. After observing data obtained 
at 40 hours, it was concluded that the data might not be valid due to cymene 
oil evaporation. Thus it was agreed upon to discard the data obtained at 40 
hours, except for using it to obtain the growth rate value for the model. 

Germination tube length was measured by including both the spore diam- 
eter and the germination tube. By observing non-germinated spore lengths 
over the course of the measurements, it was adjudicated that a non-germinated 
spore has diameter of 5 /ims or under. The data was further treated based 
on this assumption to obtain spore germination percentage for each concen- 
tration at each time point. Absolute percentage germination and percentage 
germination relative to ethanol were plotted against the log of concentra- 
tion at each time point in order to observe how percentage germination has 
changed throughout the time and how each cymene concentration has influ- 
enced spore germination compared to ethanol and to other cymene concen- 
trations. Similarly, mean absolute germ tube length and relative germ tube 
lengths were plotted against log of concentration to determine if some of the 
cymene concentrations have inhibitory or stimulatory effects on growth of 
Beauveria. This behavior would not be surprising because spores are coated 
with a strongly hydrophobic layer composed of hydrophobin that repels wa- 
ter molecules. Hence, cymene concentrations might have different impact on 
the hydrophobic behavior and the water absorption itself. Statistical analy- 
ses were performed using Minitab, Matlab, and R mathematical softwares in 
order to understand fungal growth behavior in presence of different cymene 
concentrations and acquire an idea about the possibilities of future modeling 
strategy. 

2.1. Model assumptions 

Two important assumptions were used in our model of cymene effects on 
spore germination and fungal growth. The first assumption is that there is an 
the inversely proportional relationship between the growth rate of germ tube 
and the lag time period. Snow (1949) assumed that this observation holds 
true for both a single spore and the whole population [3J. The relationship 
is demonstrated by the equation: 

rate * lag = k (2) 

where k is a constant. Dantigny (2007) verified this relationship for pop- 
ulations of eight different species of fungus [I]. Our study will test this 



5 



relationship with individual spores. 

The second assumption of our model is the dilution of water by cymene. 
Our hypothesis is that to germinate, a spore absorbs water within its vicinity 
to reach a critical volume and then germination tube starts to form. The 
time at which it reaches this critical volume is its lag time. The volume of 
spore at time t is determined by the following equation: 

^ = (1 + _H_)(l_ * )„ (3) 
dt (3 + x jy 5 + x' v ' 

In the formula, V c is the critical volume and a is the water absorption 
rate for the control. The quantity x is a theoretical ratio of volume cymene to 
volume water available in a spores vicinity. The relationship between cymene 
concentrations and ratio x is unknown. However, as cymene concentration 
decreases ten folds, ratio x also decreases ten folds. This ratio is for control 
and greater than for the other cymene concentrations. For concentration 
500 mM, ratio x = 0.7 is chosen arbitrary to yield one of the best fit. Further 
chi-square test can be used to evaluate the goodness of fit. 

A spore starts to germinate when its volume equals the critical volume, i.e. 
when y- = 1. At control, x = 0, thus, water absorption rate is a = r _1 (/i _1 ) 
where r is the lag time of an individual spore. When cymene was added, 
it had both stimulatory and inhibitory effects on water absorption rates, 
which in turn influenced lag time and percentage germination. These effects 
are apparent in the percentage germination data at 16 hours (Fig [3]). As 
cymene concentration increases from mM (control) to 5 mM, percentage 
germination increases. At 50 mM, percentage germination started to drop 
and there was no germination at 500 mM. 

Water absorption rate for an individual spore at different cymene con- 
centration is calculated by: 

a/ =(l + JZ?_)(l_^_) a (4) 
v (3 + x' K 5 + x' y ' 

where a! is water absorption rate corresponding to ratio x and a is the rate 
at x = 0. When spore reaches its critical volume at lag time rf, we have 
a/rf = 1, and thus: 

T ' ~ (1 + —) ^ 

where r is lag time at control. As water absorption rate increases, lag time 
decreases, the spore takes less time to germinate. Therefore, if cymene en- 
hances water absorption rate, it will increase percentage germination as well 
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Figure 3: Percentage germination was plotted on log cymene concentrations at 16 hours 
for control and six concentrations of cymene: 0.005 mM, 0.05 mM, 0.5 mM, 5 raM, 50 
mM and 500 mM. For control, -3 serves as log of concentration. From mM (control) to 
5 mM, cymene has stimulatory effect on percentage germination. Percentage germination 
starts to drop at 50 mM, and is at 500 mM. 

and vice versa. Water absorption rate and percentage germination vary in 
the same direction upon changes of ratio x. 

In Equation [3} parameters /3, 7, 5 are extracted from 16 hour germination 
percentage data (Fig. [3]). At concentration 5 mM, the value of percentage 
germination was highest and was approximately one and a half times the 
control. This observation was used to estimate the peak factor 7 = 0.7. 
Percentage germination at 0.5 mM is about midway between control and 
5 mM. This concentration is corresponding to x = 0.0007 and thus /3 is 
estimated to be 0.001. This value is the concentration at which 50 % of the 
population is affected by the stimulatory effect of cymene (EC 50 value). At 
50 mM, percentage germination is approximately midway between the values 
at 5 mM and 500 mM. Concentration 50 mM is the same as x value 0.07. 
The parameter 5 is estimated to be 0.1, the EC$o of inhibitory effect on spore 
germination. 
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Figure 4: Natural log of net tube lengths are plotted against time for control. Each data 
point is the mean of net tube lengths with corresponding error bars. The fitted curve has 
slope and intercept estimated by least square linear regression. 

3. Results 

3.1. Growth rate under control condition 

In control, no cymene was available. Data were obtained for 12, 16, 23 
and 40 hours as described above. From observations, spore diameter averaged 
about 5 jums, ranging from three to seven micrometers. Thus, 5 /ims was 
chosen as a threshold for spore diameter. Any measurement under 5 /ims 
was considered an ungerminated spore. These spores were not taken into 
account when growth rate was calculated. For germinated spores, 5 /ims was 
subtracted from observed length to obtain net length tube without spore 
diameter. 

Net tube lengths were exponentially proportional to time. Natural log- 
arithm transformation was performed on tube lengths. Log tube lengths 
increased linearly as time increased. Least square linear regression was used 
to estimate parameters of the model. The fitted equation is: 

log (net length tube) = 0.70186 + 0.07946 * time (6) 
net length tube = 2.02e (0 - 07946 * time) (7) 
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The estimated growth rate was 0.07946 h^ 1 with p- value 0.0299 (< 0.05). 
R square value was 0.9411, which means 94.11% of the variability of net tube 
lengths can be explained by time, i.e., the goodness of fit is substantial. Fig. 
|4| The intercept value determined the initial value for lengths of germination 
tube. The initial length value did not have any real biological meaning. At 
time 0, the initial length is 2.02 (= e 0,702 ) /mis as calculated. However, the 
model summary indicated that the intercept value 0.70186 is not statistically 
significant (p-value is 0.1855 > 0.05). Thus, the initial length values will be 
adjusted later to reach agreement between experimental and theoretical data 
for the control. 

3.2. Lag time distribution 

Lag time is defined as the time a spore spent in its inactivated stage. Per- 
centage germination was the ratio of the number of observations with lengths 
greater than 5 /ims to the total number of observations. At 4 and 8 hours, no 
germination was observed. At 40 hours, 96% of observations are greater than 
5 /ims. However, we observed that all spores have germinated at this time. 
This discrepancy happens because the spore diameters ranged from three 
to seven micrometers and five was chosen as threshold for the convenience 
of data treatment. Percentage germination versus time was plotted for the 
control case at four different time points 12, 16, 23 and 40 hours. Dantigny 
et al. (2007) assumed that fungal spore germination lag times have normal 
distribution. In this research, we also used the same assumption. Mean of 
lag time distribution was defined as the time at which half of spores have 
germinated i.e. when the percentage germination was 50%. Heuristic esti- 
mation of mean and standard deviation of lag time distribution used of the 
68-95-99.7 rule of normal distribution. Suppose variable x has normal distri- 
bution then about 68% of the values of x are within one standard deviation 
away from the mean, 95% within two standard deviations and 99.7% within 
three standard deviations. This rule is illustrated in Fig. [5j 

From plot of percentage germination versus time, the times at which the 
percentage germination were 16%, 50% and 84% were determined (Fig. |6|. 
Corresponding times are 10.15, 14.25 and 22.5 hours. Thus, estimated mean 
was 14.25 hours and estimated standard deviation was 6.125. 

Non linear regression was performed using R statistical package. The cu- 
mulative distribution function of normal distribution was the fitting function. 

W = ^(l + erf(^)) (8) 
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Figure 5: The 68-95-99.7 rule of normal distribution (adapted from Google Images) 
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Figure 6: The 68-95-99.7 rule was used to estimate mean and standard deviation of lag 
times. Horizontal lines at 0.16, 0.50 and 0.84 were drawn to estimate the times at which 
the percentage germination are 16%, 50% and 84%. 
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Figure 7: Least square nonlinear regression yielded a fitted curve to percentage germination 
plot. Cumulative distribution function is used as fitted function. 

where /i is the mean, a is the standard deviation and erf is the error function. 
Least square estimation of the mean was 14.6285 hours and standard devia- 
tion was 6.167 hours, very close to the heuristic estimations. The estimated 
mean was statistically significant. Experimental data and fitted curve for 
control were shown in Fig. [7j 

3.3. Lag time and growth rate relationship 

As mentioned in Materials and Methods, a major assumption of this study 
was the inversely proportional relationship between lag time and growth rate 
lag * rate = k. The product of the mean lag time and the estimated rate 
found above are used as the constant k. This constant does not have a unit. 

lag * rate = k = 0.07946/i * 14.6285/T 1 = 1.162381 (9) 

A thousand lag times were randomly generated from normal distribution 
with mean of 14.6285 hours and standard deviation of 6.167 hours by R. 
Some of the lag times generated from this distribution have negative values. 
This shows that normal distribution inappropriate for lag time distribution. 
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Figure 8: (a) Histogram of lag time with normal distribution with mean of 14.6285 hours 
and standard deviation of 6.167 hours. Cutoff point is 4 hours, (b) Histogram of growth 
rates calculated from lag time distribution. 



For purpose of this study, we introduced a cutoff point at 4 hours. This 
cutoff point is chosen arbitrarily to make the simulated data of percentage 
germination for control agree with experimental data. Each lag time corre- 
sponded to a spore. Growth rates were calculated from lag times using the 
relationship found in Equation [9] Distributions of lag times and growth rates 
for control are presented in Fig. [8} After growth rates were obtained, lengths 
of germination tubes were calculated using the formula: 

tube length(t) = L if t < lag (10) 
= L * e Tate{t - la9) if t > lag (11) 

Initial lengths Lq was chosen to reach agreement between lengths of ger- 
mination tubes from experiments and from model. In this study, Lq = 5/ims 
was chosen to yield that agreement. Growth curves of three randomly chosen 
spores were shown in Fig. [9] Lengths of germination tubes were calculated 
from these curves. Histograms of tube lengths from experiments and models 
were shown side by side for comparison in Fig. [TojDescriptive statistics of 
these data were presented in Table 1. 

The descriptive statistics and histograms show that the model preserves 
the positive skewness of experimental data. This means the distribution had 
long right tail and more density of data on the left. Means and standard 
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Figure 9: Three spores with three different lag times are chosen randomly from the gener- 
ated populations. Their growth rates were calculated from the relationship with lag times. 
Three growth curves were presented in the plot. 



Table 1: Descriptive statistics, comparing theoretical and experimental data for germ tube 
lengths. 





12 hours 


16 hours 


23 hours 




Expt. 


Mode 


Expt. 


Model 


Expt. 


Model 


Min 


0.833 


5 


1.786 


5 


2.036 


5 


Max 


23.324 


50.67174 


45.18 


94.16197 


63.116 


147.7996 


Mean 


4.916732 


6.414679 


11.07827 


9.389155 


13.99536 


15.83705 


Std.dev 


3.751783 


4.168555 


7.982262 


11.24463 


9.794145 


19.25989 


Median 


4.165 


5 


9.036 


5.342344 


13.234 


8.849472 


Skewness 


2.11468 


4.641366 


1.048792 


4.397724 


1.584833 


3.460279 
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Figure 10: Distribution of germ tube lengths at three experimental time periods for 
control. No treatment to data was applied. Bin size is 5 /jms. 
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Table 2: Percentage germination at 12, 16 and 23 hours for experimental and theoretical 
data. 





12h 


16h 


23h 


Experimental 


0.2926829 


0.6613546 


0.8571429 


Theoretical 


0.3179433 


0.5708290 


0.91290661 



deviations also agreed with experimental data. 

Percentage germinations for the model at 12, 16, and 23 hours were cal- 
culated. At time t, percentage germination was the ratio of the number of 
spores with lag times greater than t to the total number of spores. Percent- 
age germinations were compared between experimental and theoretical data 
in table 2. The fact that theoretical data agreed with experimental data 
verified the relationship between lag times and growth rates of the fungus B. 
bassiana 

3.4. Simulation of effect of cymene on percentage germination and growth 
of B. bassiana 

Dilution effect assumption was explained in Material and Method section. 
Model equation was constructed and parameters were estimated. 

d M = (1 + _2E_ )(1 _ » )a (12) 

alt v f3 + x JK 5 + x' v ' 

(3 = 0.001 (13) 

7 = 0.7 (14) 

5 = 0.1 (15) 

The equation and estimations of the three parameters were used to cal- 
culating percentage germinations and tube lengths at 12, 16 and 23 hours. 



Experimental and theoretical values are compared in Fig. 11 In the plot of 
tube lengths, experimental and theoretical were treated so that only germi- 
nated spores were kept. The means of net tube were plotted against loglO 
of concentrations. Theoretical data fit experimental data. Further statistical 
analysis is needed for evaluating the goodness of fit. 

Histograms produced by theoretical data preserved some important fea- 
tures of experimental data. Histograms had long right tails. This means large 
portion of population had short lengths between and 40 /ims and fewer had 
lengths over 40 /xms. As the concentration of cymene increases from to 5 
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Figure 11: Percentage germination and tube length at 12, 16 and 23 hours from experi- 
mental and theoretical data plotted side by side to show the fitness of the model. The blue 
dots are experimental data and the blue lines are theoretical. Tube length plots are shown 
on the left and percentage germination plots on the right. 
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Figure 12: Histograms of experimental and theoretical tube lengths (untreated data) plot- 
ted side by side for further comparison. Bin size is 5 /xms. Histograms are arranged in 
order of increasing cymene concentration. Heights of the bins are probabilities. 
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mM, the growth of germination tubes are stimulated by cymene. This is 
shown in the histograms by the shifting of distribution to the right. At 0.5 
and 5 mM, the highest bin was 5-10 /mi instead of 0-5 /mi as in preceding 
histograms. The widths of histograms decreased at 50 mM. Compared to 
5 mM, the width was slightly reduced. At 500 mM, there are only lengths 
between and 5 /mis, corresponding to ungerminated spores. 

4. Discussions 

Smaller concentrations of cymene have a stimulatory effect on growth of 
Beauveria bassiana, while greater concentrations have an inhibitory effect. 
Similarly, smaller concentrations of cymene attract insects, while greater 
concentrations of cymene act as insect repellant jH]- Beauveria bassiana is 
an effective entomopathogen, the concentrations that attract the insect may 
induce germination of the spores, thus allowing the fungus to locate a suitable 
host. Cymene is volatile and evaporates in a couple of hours, and thus a seed 
is protected at initial stage until B. bassiana spores germinate. Germ tubes 
in our experiment grew in all directions randomly, searching for nutrients 
that were not provided. If there are no nutrients, the fungus will not survive. 
Since plants provide food for fungi, fungi have to reach a plant as soon as 
possible in order to survive. Thus, appropriate cymene concentration should 
be present that will allow fungi to inhibit the plant and that will also repel, 
rather than attract insects. 

Our study used the assumption that lag times have normal distribution 
that was demonstrated by Dantigny (2007). This distribution produced neg- 
ative lag times in the simulation. To avoid this, Dantigny (2007) introduced 
skewness into normal distribution. Another way is to find a better fitting 
distribution, ones with all positive values; this will be one of our future di- 
rections. The percentage germination over time were used in both Dantigny's 
paper and this study. While Dantigny used Gomperzt and logistic equation 
to fit the germination data, we used the cumulative distribution function of 
normal distribution. 

A major difference in the approach used by Dantigny (2007) and the one 
in this study is that he measured the radial growth of fungal colonies. These 
data were used to calculate the growth rate of the whole population. This 
was not possible in our study because B. bassiana produced multiple colonies 
on petri dishes, making it impossible to measure radial growth. We chose 
therefore to measure germ tube lengths. 
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A primary assumption for Dantigny's and our study used was that there is 
an inverse proportionality between growth rate and lag time. While Dantigny 
tested this assumption on the whole fungal population, our study verified this 
relationship for the germination of single spores. 

This research is important because it allows a relatively simple mathe- 
matical model of the development of spore germ tubes in the absence and 
presence of a natural biopesticide. This model may prove applicable to eval- 
uation of biopesticides. 
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